UWThPh-2003-42 



Leptogenesis in seesaw models with a 
twofold-degenerate neutrino Dirac mass matrix 

Walter Grimus* 

Institut fur Theoretische Physik, Universitat Wien 
Boltzmanngasse 5, A-1090 Wien, Austria 

Luis Lavoura** 

Universidade Tecnica de Lisboa 
Centro de Ffsica das InteracgSes Fundamentais 
Instituto Superior Tecnico, P-1049-001 Lisboa, Portugal 

10 May 2004 

Abstract 

We study leptogenesis in two seesaw models where maximal atmospheric neu- 
trino mixing and U e 3 = result from symmetries. Salient features of those models 
are the existence of three Higgs doublets and a twofold degeneracy of the neutrino 
Dirac mass matrix. We find that in those models both leptogenesis and neutrinoless 
double beta decay depend on the same unique Majorana phase. Leptogenesis can 
produce a baryon asymmetry of the universe of the right size provided the mass 
of the heavy neutrino whose decays generate the lepton asymmetry is in the range 
10 11 -10 12 GeV. Moreover, in these models, leptogenesis precludes an inverted neu- 
trino mass spectrum since it requires the mass of the lightest neutrino to be in the 
range 10 _3 -10 -2 eV. 
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1 Introduction 



Experimental cosmology has witnessed spectacular progress during the last few years. In 
particular, the WMAP experiment PQ has determined with fantastic precision the baryon 
asymmetry of the universe, which is given by the ratio of baryon number over the number 
of photons, experimentally measured to be 

Vb = = 6.5i& x 10- 10 , (1) 

where n# is the (present) baryon density of the universe, n§ is the anti-baryon density 
and n 7 is the (present) density of photons. This value of tjb is in excellent agreement with 
the one inferred from big bang nucleosynthesis |2J. 

Another field which has lately witnessed outstanding experimental progress is neutrino 
masses and lepton mixing. On the one hand, the first results of the KamLAND experiment 
[3] have conclusively proved that solar neutrinos oscillate; a global analysis of all solar 
neutrino results, including the recent SNO measurement [I] and also the KamLAND 
result, gave a mass-squared difference [I] 

Am 2 Q = m\-m\ = 7.1±5;| x 10" 5 eV 2 (2) 

and a large but non-maximal mixing angle 

9 = 32.5+ 2 ^ degrees, (3) 

where the errors reflect la constraints in the two-dimensional 0— Am^ region. Note that 
tan# is the ratio, in the decomposition of the electron neutrino v e , of the amplitude for 
the heavier neutrino V2 over the amplitude for the lighter neutrino v\\ tan# = \U e 2 /U& \, 
where U is the lepton mixing matrix. On the other hand, the Super-Kamiokande exper- 
iment [5] has shown that atmospheric neutrinos oscillate with a mass-squared difference 
in the 90% CL range 



1.3 x 10" 3 eV 2 < Am 



atm 



m 3 — m-. 



<3.0xl0^eV 2 , (4) 



with best-fit value Am 2 tm = 2.0 x 10 3 eV 2 , and a (most likely) maximal mixing angle: 

2 2fl atm ee 4 {U^l 2 (l - p^l 2 ) > 0.90 (5) 



sin 



at 90% CL, the best-fit value being exactly 1. Finally, the CHOOZ experiment [E] and all 
other neutrino oscillation data yield the upper bound \U e 3\ 2 ^ 0.054 at 3cr [7j. For recent 
reviews on neutrino oscillations see jSj. 

These experimental developments invite a renewed interest of theorists for leptogenesis 
[HI QS]. This is the possibility that the baryon asymmetry of the universe has been 
generated through the standard-model sphaleron transmutation of a previously existing 
lepton asymmetry, which in turn was generated at the decay of the heavy neutrinos 
involved in the seesaw mechanism jllj . In the standard version of that mechanism one 
introduces three gauge-singlet right-handed neutrinos v R . Let us define 



v' L ee Cul, (6) 
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where C is the charge-conjugation matrix. Then, the mass terms for the neutrinos are 
given by J2] 

CumasB = -VrM D Vl ~ V L M\jVr -^VrMrCv^ + ^V^C' 1 M r Vr (7) 

where Mr is a symmetric matrix. If the eigenvalues of MrM^ are all much larger than 
the eigenvalues of MdM],, then the approximate Majorana mass matrix for the light 
neutrinos is given by 

M v = -M T D M^M Dl (9) 

while the Majorana mass matrix for the heavy neutrinos is approximately equal to Mr 
[T5] . In the weak basis where the mass matrix M^ of the charged leptons is diagonal, 
Mi = diag (m e , m M , m T ), one has 

U T M V U = diag (m l5 m 2 , m 3 ) , (10) 

where m lj2 3 are real non-negative and C/ is once again the lepton mixing matrix. 

Since leptogenesis needs CP violation, an intriguing question is whether there is a 
connection between the CP violation at the seesaw scale and the one at low energies. In 
the most general case, the answer to this question is negative [TH ll5|. However, it is easy 
to find scenarios where such a connection exists — see, for instance, [THl I17j; in minimal 
scenarios only two heavy Majorana neutrinos are required ^Sl El 1201 E] Inspired by 
grand unified theories, it is also quite common in studies of leptogenesis to assume hierar- 
chies in the neutrino sector — see [221 ESI 12H I2H] and references therein. In particular, one 
may assume that the Dirac neutrino mass matrix Mo is strongly hierarchical, i.e. that 
the eigenvalues |a| 2 , \b\ 2 , |c| 2 of MjjMjj satisfy \a\ <C \b\ <C |c|, and subsequently recon- 
struct the masses of the heavy Majorana neutrinos from the low energy data [221 123] ■ The 
assumption of a hierarchy in Mq is justified by the relationship, existing in some grand 
unified theories, between Mp and the up-type-quark mass matrix, and by the fact that 
the latter matrix is known to be strongly hierarchical. 

In this paper we take a different departure and start from the fact that atmospheric 
neutrino mixing is maximal (or nearly maximal) , which suggests an alternative possibility 
[2H1 [23 I2H] • Since experimentally ([/^l ~ \U T s\, there may exist in nature a fi-r inter- 
change symmetry. We know that ^ m T , hence the /i-r interchange symmetry must be 
broken in the charged-lepton sector, but it may be kept intact in the neutrino sector. This 
can be achieved through the introduction of three Higgs doublets, one of them ((pi) with 
Yukawa couplings to the neutrino singlets and to the charged-lepton singlet e#, and the 
other two (0 2 an d cj) 3 ) with Yukawa couplings only to the charged-lepton singlets fj,R and 
Tr. Under the interchange fi <-> r the doublet </> 2 remains invariant while 03 changes sign; 
this leads to m M 7^ m T . On the other hand, the neutrino Dirac mass matrix is twofold 
degenerate because 4>i is invariant under <-> r: 

Md = diag (a, 6, b) . (11) 



3 



A crucial feature of these models is the existence of some other symmetry — either the 
continuous lepton-number symmetries |2EJ or a discrete symmetry |2Zj — which forces 
and the charged-lepton mass matrix Mi to be simultaneously diagonal. These other 
symmetries are allowed to be softly broken, hence the right-handed-neutrino Majorana 
mass matrix is non-diagonal and has the form 





( 


in 


n 


n 


M R = 




n 


V 


q 




\ 


n 


Q 


V 



due to the /i-r interchange symmetry. It is this matrix Mr which produces lepton mixing. 

The neutrino sectors of the model of j2H| and of the D4 model of |2Z| are both 
characterized by equations (jll)) and (fT2*|) : the -D4 model is more constrained than the 7L>2 
model since it has q = in Mr. Note that the CP-violating phase analogous to the CKM 
phase is absent from the models under discussion, because U e 3 = 0. This follows easily 
from equations (fTTj) and (fT^j) . since 

/ x y y \ 

M u = \ y z w \ (13) 
\y w z J 

in the basis in which Mi is diagonal. Thus the only sources of CP violation in the leptonic 
sector are the two physical Majorana phases in U . 

The purpose of this paper consists in analyzing leptogenesis in the models of [2*d1 12*7] . 
In particular, we will show that they have the following properties: 

1. Leptogenesis is a viable scenario. 

2. Correctly reproducing rjB constrains the spectra of both the light and the heavy 
neutrinos. 

3. Only one of the two Majorana phases is responsible for leptogenesis, and that phase 
is also the only one which appears in the effective mass | (ttt.) | for neutrinoless (3/3 
decay. 

In section |21 we review the computation of t]b from the knowledge of Mjj and Mr, with 
emphasis on the three-Higgs-doublet structure of our models. We proceed in section El to 
derive the relevant analytic formulae for the diagonalization of Mr and Ai u , in order to 
calculate t]b- We apply those formulae in section HI to study the variation of tjb with the 
parameters of the models. In section El we draw our conclusions. An appendix contains 
calculational details related to section El 

2 Baryogenesis from leptogenesis 

The "natural" basis for our models is given by diagonal matrices Mb and Mi while Mr is 
non-diagonal — see equations (|llj) and ()12|). However, the basis in which the leptogenesis 
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formalism is established is the one where Mi and the mass matrix of the right-handed 
neutrino singlets are diagonal; the latter matrix is then 



M R = diag (M ls M 2) M 3 ), 
with real non-negative diagonal elements. Defining a unitary matrix V by 



we find, using equation (JJJ), 



V T M R V = M R , 



M' D = V T M D 



(14) 



(15) 



(16) 



for the neutrino Dirac mass matrix in the leptogenesis basis. For the actual calculation 
of t]b one needs 

(17) 



R = M'nM'J = V T M D MlV 



We assume, for the masses of the heavy Majorana neutrinos iV^^, that Mi <C M%£. 
Then, the CP asymmetry ei produced in the decay of N 1 (the heavy neutrino with mass 
Mi) is HDIESI 



1 



?7T It'll R\ 



M 



1 j=2 



M{ 



Im 



[R 



where V\ denotes the vacuum expectation value (VEV) of 4>\. The function / is given by 



f(t) = Vt 



2-t N1 t 

1 + 1) In 



For t ^> 1, we have 

/(*) 



'1/2 



t 



-3/2 



6 



13 
12 



-5/2 



l+t 



19 
20 



-7/2 



(19) 



(20) 



Thus, / (£) is negative. 

The leptonic asymmetry produced through the decay of Ni is written as [TUl E01 IUj 



n L - n L 



(21) 



where is the lepton density, is the anti-lepton density, s is the entropy density, K\ 
is the dilution factor for the CP asymmetry e\ and g*i is the effective number of degrees 
of freedom at the temperature T = Mi. The effective number of degrees of freedom is 
given by (see for instance 



7 

9* = 9j + g 9k- 

j'=boson fc=fermion 



[22) 



In the SU(2) x Z7(l) gauge theory with three Higgs doublets and supplemented by the 
seesaw mechanism one has 



9*i 



28 



x 90 



SM 

5 



+ 8 + - x 2 = 116.5, 



(23) 



where the terms within the brackets are the Standard Model contributions and the last 
two terms in the sum take into account the two additional Higgs doublets and the lightest 
heavy neutrino Ni, respectively. 

The baryon asymmetry Y B produced through the sphaleron transmutation of Yl, while 
the quantum number B — L remains conserved, is given by 

v 00 v 8iV> + 4N H 

Yb = u with to = , 24 

lu-1 22N F + 13N H ' v ; 

where Np = 3 is the number of fermion families and Nh is the number of Higgs doublets. 
This relation derives from the thermal equilibrium of sphalerons for 10 2 GeV^T^lO 12 GeV 
[TUJ [BTJ . Note that we must impose the condition Mi < 10 12 GeV, otherwise Y L could be 
erased before it transmutes into Yb- From equation (J24j) . uj = 12/35 in three-Higgs- 
doublet models. 

Next we discuss the relation between Yb and rjs, where the latter quantity is the 
ratio baryon-number density over photon density. Note that this is the present ratio, and 
rjB has last changed at the time of e + e~ annihilation, at which time the photon density 
(temperature) has increased relative to the neutrino density (temperature). On the other 
hand, Yb did not change since its generation — the baryon number per comoving volume 
and the entropy per comoving volume remain constant. Thus we have 



s 

Vb= — 

n 7 



Y B , (25) 



o 



where n 7 is the photon density. The index denotes present time. We know that [32] 

s = ^g* T 3 and n 7 = -|c(3)T 3 , (26) 

where T is the photon temperature and ( (3) ~ 1.20206, ( being the zeta function. 
At present only photons and light neutrinos are relevant for s, since all other particles 
have annihilated apart from tiny remnants — actually, we are just discussing the baryon 
remnant. Thus 

7 4 43 

<fc = 2 + -x6x— = -, (27) 

where we have taken into account that the neutrino temperature after e + e~ annihilation 
is a factor of (4/H) 1 ^ 3 lower than the photon temperature T. One thus obtains [201 HH] 

7T^ 43 

7.0394 and r] B = 7.0394 Y B ~ -3.15 x 10" 2 K X e x . (28) 



C(3) 495 



We have used the values for three-Higgs-doublet models. 

We now turn to the dilution factor, which is approximately given by |30| EE1 



where 

Ki = (30) 
ill 
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In this equation, I\ is the decay width of Ni, given at tree level by 



r. = (3D 

07T \Vi\ 



and Hi is the Hubble constant at temperature T = M±, 



M 2 

H 1 = 1.66Jg^—^, (32) 

J W Planck 



where M P i anck = 1.221 x 10 19 GeV. Thus 



K X = Mp ' anckjRll 2 . (33) 

lM+Zg^&K | vi\ 2 Mi 



Equation (JHJ) holds for 10 < K x < 10 6 [HI E3 Numerically one obtains 

895.6 i? n Zl74GeV\ 2 , . 

^w^n^H ■ (34) 

For g*i we have used the value of equation (|23[). 

With equations (fTHJ) . (J25|) . ()29|) and ()33|) one obtains, after a numerical evaluation, 

^--1.39x10 9 ^/ _L L — . (35) 



(InK.f 5 p 2 J \Ml) (i? n ) 2 10 11 GeV 

This equation, which we shall use in conjunction with equation ()34j) . clearly indicates 
the desired order of magnitude of M\. Equation (|33|) shows that tjb is approximately 
independent of \v\\. Only the logarithm of K\ in the denominator makes t]b dependent 
on \v\\. Also notice that the minus sign in equation ()35|) cancels with the negative sign 

of/(0- 



3 The baryon asymmetry in our models 

It remains to calculate the quantities Im [(i?ij) 2 (j = 2,3) and Rn in our models. We 
will see that the Z2 model allows a full analytical calculation of these quantities and, 
according to the discussion at the end of section ^ these quantities in the D4 model are 
special cases of those in the 1^2 model. 

We remind the reader that, in our models, the effective light-neutrino Majorana mass 
matrix Ai u is as in equation (|13|). while the right-handed-neutrino Majorana mass matrix 
Mr has the same form as hA v and is in equation ([12ft . The matrix M. v is diagonalized by 
U, see equation (fTOjl. while M R is diagonalized by V, see equations (fT4*j) and (fT5|). The 
analogies between M. v and M R and between U and V are used repeatedly throughout 
this paper. 
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It can be shown that, because M. v is of the specific form in equation (JT3J), i.e. with 
(•M„) = (Ai v ) TT and (M. u ) = (A4 u ) eT , the matrix U can be parametrized as 



U 



diag (l, 



/ -c s 

s/V2 c/y/2 l/y/2 diag(< 
V s/y/2 c/V2 -l/y/2 



(36) 



with c = cos# and s = sin^. 1 We assume, without loss of generality, that 9 belongs to 
the first quadrant while m 2 > m\. The phase a is unphysical; the only physical phases 
are the differences 

A = 2(0! -ft) (37) 

and 2 (ft — ft). The matrix .M;, has six physical parameters: the moduli of x, y, z and w 
and the phases of zw* and y 2 x*z*. 2 It is better to use as physical parameters the moduli 
of z + w and z — w instead of the moduli of z and w, and the phases of (z — w) (z + w)* 
and y 2 x* (z + w)* instead of the phases of zw* and y 2 x*z*. The six parameters \x\, 
\y\, \z + w\, \z — w\, arg [y 2 x* (z + w)*] and arg [(z — w) (z + w)*] correspond to the six 
observables mi, m 2 , m 3 , 9, A and 2(ft — ft). 

Since the matrices Mr and A4 U have the same form, the matrix V which diagonalizes 
Mr has the same form as the matrix U which diagonalizes M. v \ 



V 



diag (l 



, e lx , e lx 



t -d s' 

s'/V2 d/y/2 l/y/2 ] diag(e i7l ,e i72 ,e 
V s'/y/2 d/y/2 -l/y/2 



'73 



(3? 



with d = cos 9' and s' = sin^'. We assume, without loss of generality, that 9' belongs 
to the first quadrant while M 2 > M±. Like in the previous paragraph, there are six 
observables originating in Mr: Mi, M 2 , M 3 , 9', 2(71 — 72) and 2(71 — 73). 

We proceed to the calculation of R — see equation (fTTj) . Using equations (fTTj) and (JHBJ), 
we obtain 



R 



V T M D M ] D V* 



( 



ds' 



|a| 2 c' 2 +|6| 2 S ' 2 



lal 2 ) e i(72 " 7l) 



z's 1 (|6| : 



2> \ p»(7i-72) 



I 1 2 /2 , 1 7 1 2 /2 

|a| s + \b\ c 



\ 


\bf ) 



(39) 



The fact that 



R 



13 



2:-! 







(40) 



implies that in this model the third heavy neutrino has no bearing on leptogenesis — even 
if its mass M 3 happens to be lower than M\ and M 2 , the masses of the other two heavy 
neutrinos. 

The matrices M. v and Mr are related through equation with Mb given by equa- 
tion pijl . That relation contains only two extra parameters: \a\ and |6|. This means that, 

1 Equation (|36|l is clearly not the most general parametrization for a unitary matrix, rather it is a 
consequence of the specific form of M. v in equation (|13|) . 

2 The physical phases in M. v are the ones of its Jarlskog invariants. For instance, arg (zw*) is the phase 
of [M v )„ [M v )* {M v )* and arg (y*x*z*) is the phase of (M v ) eu [M v ) {M v )le {M v )* 
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out of the six observables originating in Mr, only two can be considered as independent of 
the six observables originating in M. v . We shall select M 12 , together with the observables 
originating in Ai u , as the basic observables of the theory, and express the four remain- 
ing ones in terms of these. It will turn out that 2 (fli — f3^) and play no role in the 
computation of t]b', the basic observables that we need for that computation are 



mi 2, Mi 2, and A. 



(41) 



With the aim of using the experimental information on AttIq and Am 2 tm , we express 



m 2 



mi + Am® and = y mf + a Am 



) 2 

""atm 



<7 



±1) 



(42) 



as functions of m\. For m% there is a twofold ambiguity stemming from the two possible 
neutrino spectra: normal spectrum mi < m 2 < m 3 (a = +1) and inverted spectrum 

m 3 < mi < m 2 (c = — 1)- The inverted spectrum is only allowed for mi > y'Am^; the 
experimental result for the atmospheric mass-squared difference — see equation (J3J) — then 
requires mi > 0.04 eV, which is rather large. 

As shown in the Appendix, the parameters \x\, \y\ and \z + w\ are given in terms of 
the observables in equation (}4Tj) by 



Fl 

\y\ 

\z + w\ 



2 1 2 iA 

c nil + s m 2 e 



cs 
V2 



nil — 1^2^ 



2 1 2 iA 

s nil + c m 2 e 



We may then compute 
B 
C 



4mim2MiM 2 \y\ 



2 2 
m 1 m 2 



(M x 2 + M 2 2 ) 



m\m\M\Ml \x\ \z + w\ 



These are the coefficients of the equation 



|x| 4 |6| 8 + fi|a;| 2 |fe| 4 + C = 0, 



(43) 
(44) 

(45) 

(46) 
(47) 

(48) 



which is derived in the Appendix and allows one to compute the parameter \b\ as a function 
of the observables in equation (}4~Tj) . Indeed, 



-B±VB 2 -4:C 
2 Id 2 



Moreover, the parameter \a\ is given by — see equation ()A32|) 



-B T VB 2 - AC 
2\z + w\ 2 



(49) 



(50) 



and the mixing angle 9' is given, as a function of the observables in equation (|41|). by — see 
equation (jA34j) 

j 2 



l2 

C — S 



±VB 2 - AC 
m\m\ (Mf - Mf) 
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(51) 



Thus, using m 12 , M 12 , 6 and A as input, we are able to compute, with a twofold ambigu- 
ity, |6|, \a\ and 9'. The twofold ambiguity corresponds to the interchanges \m\ <-> \p + q\ 



and c' 2 as is made clear in the Appendix. 

One must check that the condition B 2 — AC > is respected. This condition leads to 
a lower bound on mi; in the limit M 2 3> Mi one finds approximately 



,/2 



m x > M x 
m 2 ~ M 2 



^ -pp- sin 2 29 or mi > J AtHq — ^ sin 2 20. 



M, 



(52) 



On the other hand, it turns out that, in general, with m 2 given by the first equation 
|a| 2 and |6| 2 increase with increasing m 1 . In our models, a and 6 arise from the Yukawa 
Lagrangian 



Y 



*5, 



+ b 



V, 



VtR 



h.c. 



, (53) 



where • • ■ represents the Yukawa couplings of the right-handed charged-lepton fields. If 
we require that the Yukawa coupling constants a/v\ and b/v± should at most be of order 
1, this implies a loose upper bound on mi (when m\ — m\ is kept fixed), since \a\ and |6| 
increase with mi. For M 2 3> Mi and using the solution \b\ 2 > \a\ 2 , i.e. the upper signs in 
equations and (fSTJj). we find the approximate expressions 



miMi 1 - sin 2 (A/2) sin 2 29 
miM 2 [l - sin 2 (A/2) sin 2 29 



1/2 
-1/2 



(54) 
(55) 



where we have used m 2 — mi, valid for mi ^> y Am®. These equations may be used 
to compute the approximate upper bound on mi. We see that \b\ 2 ~ m\M 2 . Therefore, 



requiring 



1 2 5; 1 leads to 



mi 



< 



Pi I 
M 9 



(56) 



For instance, for \v\ \ = 10 GeV and M 2 = 10 13 GeV, equation yields mi < 0.01 eV. If 
we choose \vi\ = 50 GeV and M 2 = 2.5 x 10 12 GeV, we have mi < 1 eV. 

An expression for Im (-Ri 2 ) 2 in terms of our basic observables is derived in the Ap- 



pendix, see equation (jA27j) . Using equation the main result of this section is 

1 „/M 2 \/,. l9 , l9X 2 Mf M 2 (m| — mf) c 2 s 2 sin A 



rjB - -1.39 x 10 



-9 



(hi 



,3/5 



'(3)01' 



mim 2 (Mf - Mf) (i?n) 2 (10 11 GeV) 



Here, a convenient expression for R n — see equation (}3T?j) — is 



|a| 2 + |6| 2 + (|a| 2 -|6| 2 ) (V 



2 /2 
S 



(57) 
(58) 



where c' 2 — s' 2 is given by equation (foTj) . 
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Since / (Mf /Mf) is negative, we find from equation (J57|) that sin A > 0. Notice the 
important point that the violation of CP responsible for the generation of a non-zero t]b 
all comes from the Majorana phase A. This is the same Majorana phase entering the 
matrix element for neutrinoless double beta decay, 

IWI = l(M,)J = |s|, (59) 

which is given by equation ()43|). Thus, in these models neutrinoless double beta decay 
and leptogenesis depend on the same Majorana phase A. 

To conclude this section, we discuss the extra condition on Mr in the D4 model of 
|2"7| . As mentioned before, in that model (Mr) = q — 0. Since Mr, is diagonal, this 
leads to {M' 1 )^ = 0, or [27] 

il e 2 ^ + e *fc _ J_ e m = (60 ) 
mi m 2 m 3 

The Majorana phase 2(/3i — ^3) is irrelevant for our purposes, therefore the important 
constraint is 



S 2 ,A C 2 

— e H 

mi m2 



(61) 



1 

m 3 

This forces m 3 to be larger than both m x and m 2 (normal spectrum). The Majorana 
phase A becomes a function of mi through 



(mim 2 /ffl 3 ) -c 4 m 2 -s 4 m^ 
cosA = — ^ 62 

and through equations ()42j) with cr = +1. Thus, the -D4 model of [2Zj has one degree of 
freedom less than the Z 2 model of |2"6l . 



4 Numerical results 

In this section we shall always use 

= 33°, (63) 
Am^ = 7.1xlO~ 5 eV 2 , (64) 

the best-fit values of j3j. Then the observables in our set (|41jl which are still free to be 
chosen are mi, M± t 2 and A; we furthermore have to choose \v±\ — see equations (jMj) and 
(foTj) . The mass m 2 is fixed via equations (Jj2j) and (jHU)- 

From this input one obtains \x\, \y\ and \z + w\ using equations (j43|l - (j45|l . One then 
computes \a\, \b\ and 9', with a twofold ambiguity, from equations (f47Jj) - (f5T]) . Thereafter, 
Rn is found in equation (|5£J) and K\ is given by equation (jHl|). 

As for the twofold ambiguity in the computation of \a\ and |6|, we use the upper signs 
in equations (j4T?|) - (joT|) . Numerically, the choice of the lower signs yields a smaller t]b- 

The VEV \vi\ must be smaller than 174 GeV, in order that there is also room for 
the other two VEVs: £f=i |^| 2 = (174 GeV) 2 . Since in the models of gEl EJJ both the 
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0.001 0.01 0.1 

m, (eV) 



Figure 1: t\b as function of m\ for three different values of M\. In producing this figure 
we have chosen M 2 /Mi = 10, \v\\ = 50 GeV and A = 90°. The lowest allowed mi for this 
value M2/M1 = 10 is mi = 0.71 x 10~ 3 eV. The horizontal lines indicate the experimental 
result for t]b, as given in equation (JTJ). 

neutrino masses and the electron mass originate in Yukawa couplings to 0i, while the /x 
and t masses originate in Yukawa couplings to 02 and 03, the smallness of the electron 
and neutrino masses suggests that V\ should be relatively small. 

In figure ^ we have plotted t\b versus mi for three different values of M\. We read off 
from that figure that M\ must be larger than 10 11 GeV in order to reproduce the experi- 
mental value of t]b in equation (JTJ). Since a successful leptogenesis requires Mi < 10 12 GeV 
[TIH EH], the order of magnitude of M\ becomes quite constrained. We furthermore see 
that small values of m\ are preferred; for large values of mi, t]B becomes too small. Ac- 
tually from figure^we gather that mi cannot exceed 0.02 eV. This rules out the inverted 
neutrino spectrum if we want to accommodate leptogenesis in our model. As a function of 
mi, the maximum of the theoretical expression ([57]) for t]b is roughly at mi = 3 x 10~ 3 eV. 
If we consider t)b as a function of the CP-violating Majorana phase A, a numerical study 
shows that the maximum of T]b is attained for A close to 100°. As a function of \v\\, t\b 
increases by less than a factor of two when that VEV goes from 10 GeV to 100 GeV. 

In order to understand the dependence of t\b of equation (jBTj) on Mi j2 , it is useful to 
perform a scale transformation —> AMi^ where A is an arbitrary positive number. 
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Figure 2: t]b as a function of mi for Mi = 2.5 x 10 11 GeV and three different values of 
M 2 . We have used \v x \ = 50 GeV and A = 90°. 

From equations and (|30)1 we see that |a| 2 and |6| 2 scale with one power of A; the same 
holds for R12 and Rn. It follows that t)b also scales with one power of A; in other words, 
rjB is a homogeneous function of order one in Mi i2 . 

In figure|21we have plotted tjb versus mi for M\ = 2.5 x 10 11 GeV and different values of 
M 2 . This figure shows that fixing Mi and varying M 2 with M 2 3> Mi does not drastically 
alter t]b, except for very small values of mi. The lower bound on mi as a function of the 
ratio M 2 /Mi is also illustrated in figure El 

5 Conclusions 

In this paper we have computed the baryon asymmetry of the universe in the Z 2 model 
of |2E1 and in the -D4 model of j2ZJ- These models are characterized by a neutrino Dirac 
mass matrix Mjj with two degenerate eigenvalues, and by an interchange (Z 2 ) symmetry 
between the fi and r families. Both models predict maximal atmospheric neutrino mixing 
and U e 3 = as a consequence of their symmetries. 

We have shown that these models can easily accommodate baryogenesis via leptoge- 
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Figure 3: The effective Majorana mass in neutrinoless (3(3 decay as a function of m\. 
Starting at the upper curve and descending to the lowest of the four curves, the values 
used for A are 0°, 90°, 135° and 180°, respectively. 

nesis and reproduce the experimental value of tjb, 3 provided the mass of Ni, the light- 
est heavy neutrino whose decays generate the relevant CP asymmetry, is in the range 
10 11 GeV < Mi < 10 12 GeV. Furthermore, the mass mi of the lightest neutrino has to 
be in the range 10 -3 eV ^ mi < 10~ 2 eV. Thus, an inverted neutrino mass spectrum is 
practically excluded. 

In figure El we have plotted the effective mass probed in neutrinoless (3(3 decay. 

We see that, if our models are to have successful leptogenesis, then | (m) | is at most 0.01 eV; 
if the evidence for neutrinoless (3(3 decay of |37], with \(m)\ > 0.1 eV is confirmed, then 
leptogenesis is not enough to generate an t]b of the observed size. 

We stress that the neutrino mass matrices of our models allow an analytical calculation 
of the CP asymmetry ei of equation (fTSj) . The results of this paper for the and D4 
models may be valid in a wider context of general models with a twofold degenerate Dirac 
mass matrix M^; the reason is that the mass matrix of the light neutrinos must have the 
form in equation (|13|) if one assumes maximal atmospheric neutrino mixing and U e 3 = 0, 

3 On the other hand, as shown in |2$j . a mass hierarchy in Md requires finetuning of the masses in the 
heavy neutrino sector in order to reproduce t\b- 
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assumptions that, as experiment shows, cannot be far from true. 
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A Algebraic details 

The matrices Mr and M. v Those matrices are given in equations (|12|) and (|13)). and 
their relationship to each other is expressed by equations (0) and (JTTJ). We define 

d = 2n 2 - m (p + q) , (Al) 
/ = 2y 2 -x(z + w) J (A2) 

so that 

det M R = (q - p) d and det M u = (w - z) f. (A3) 
Then, by explicitly inverting Mr and by using equation Q, we find 

x = a 2 ^-, (A4) 
a 

— 71 

y = ab -j-' ( A5 ) 

771 

z + w = b 2 — (A6) 

and 

(z -w)(q-p) = b 2 . (A7) 
It is also useful to invert equation Q: 



Mr = -M D M~ Y M L D . (AS 



From this relation we compute 



o z + w ... 

m = a 2 — — , (A9) 

n = aby-, (A10) 

P + q = b 2 j. (All) 



In addition, we obtain equation (jA7|) again. 

The parameters x, y, z and w From equations (fTUj) . f)13|) and (J3*6*j) one may write 

/ x ye ia ye ia \ I -c s \ / -c s/v 7 ^ s/v 7 ^ \ 

ye ia ^e 2ia we 2M = s/v^ c/v^ 1/^2 /} s c/v^ c/v 7 ^ 
V ye M we 2M ze 2ia / \ c/v^ -1/^2 / \ 1/^2 -1/v^ / 

(A12) 

where 

fx = diag (m ie - 2t/31 , m 2 e- 2 ^ 2 , m 3 e- 2i/33 ) . (A13) 
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From this we obtain equations giving the parameters of M. v as functions of the physical 
observables. We find equations (pf^ - Po")) and also 



y 2 x* (z + w) 



c 2 s 2 



(mi - m 2 e iA ) 2 (c 2 mi + s 2 m 2 e~ iA ) (s 2 mi + c 2 m 2 e~ iA ) , (A14) 

(A15) 
(A16) 



2 

\z — w\ = m 3 
[z — w) [z + w)* = m 3 
In analogy to equation ()A15j) we also have, for Mr instead of M. v 

\p - q\ = M 3 . 



s 2 mi e 2i ^-M + c 2 m 2 e 2 ^-^ 



(A17) 



The parameters a and b We now express \a\ and |6| as functions of the neutrino 
masses. Using equations (|A7jl . (jA15|) and (|A17j) we obtain 

\b\ 2 = m 3 M 3 . (A18) 

2 4 

Comparing det A^ with det Mr we readily find that \a\ \b\ = m\m 2 m 3 M\M 2 M 3 . There- 
fore, with equation (|A18|) we conclude that 

l2 mvmoMxMo 

. (A19) 



m 3 M 3 

Moreover, equations (|A3|) . ()A15|) and ()A17|) lead to 

|d| = MiM 2 , 
|/| = m x m 2 . 

The imaginary part of (i?i 2 ) 2 From equation (|A14|) one derives 

2 * / i \*1 2 2 ^2 2\ ■ \ 

y x (z + w) = c s mim 2 ( m 2 — m x I sin A. 



21m 



(A20) 
(A21) 



(A22) 



Using the analogy between the matrices Ai u and Mr, and between U and V, one sees 
that 



21m 



n m* (p + q)* 



■ d 2 s' 2 M x M 2 [M\ - Ml) sin (2 7i - 2 72 ). 
The matrix R is given in equation (J39j) . Using equation ()A23|) one sees that 



Im 



(R 



12 



2 , 2 \ 2 2 Im [n 2 m* {p + g)*] 



MiAf 2 (M| - Mj 2 ) ' 
Now, from equations (|A4|) - (|A6[) one finds the relation 



2 */ , \ * | |4|,|4^*(P + ?)" 

j/ s (z + w) = \a\ | o | - 



(A23) 
(A24) 

(A25) 
(A26) 

mi m 2 (Mf - M 2 ) " " (A27) 

where we have used equations (|A18J) - (|A20J) in order to go from equation f)A2fi|) to equa- 
tion (^27) . 



Therefore, equations ()A22|) and (|A24|) give 



Im 



(\b\ 2 - 

(\b\ 2 ~ 



2 \2 m x m 2 (m\-mf) \df 2 2 . 

a ; — s ^r- — i — jcs smA 

1 ^ M!M 2 (M| - Mf) |a| 4 |6| 4 

2^2 M X M 2 [m\ -m\) 2 2 . 

c s smA, 
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The quadratic equation giving \a\ and \b\ Equations (|A18|) and (jA19|) contain m 3 
and M 3 . We shall now derive expressions for \a\ and |6| which are functions solely of the 
observables in equation (jUJ. With equations it is easy to check that 

ml + ml = \x\ 2 + 4:\y\ 2 + \z + w\ 2 . (A28) 

For M R instead of A4 U , the analogous relation is 

M 2 + Ml = \m\ 2 + 4 \n\ 2 + \p + q\ 2 . (A29) 

Using equations (lA"9"l)- (IAlll) and (IA21D . we find 

i\m\ (Ml + Ml) = \af \z + w\ 2 + 4 \a\ 2 \b\ 2 \y\ 2 + |6| 4 \x\ 2 . (A30) 



Multiplying this equation by \x\ 2 \b\ 4 and using \a\ 2 \b\ 2 = ■m\miM\Mi we finally obtain 
equations (J4*t3j) - (j4*5j) . The solutions to equation (|38|) are 



m Xb> + 9 |' = W'|t|'= - B± f^ . (A31) 

where we have used equations (jAll|) and (jA21|) . Since C = (j&| 4 \x\ 2 ) (ja| 4 \z + w\ 2 ), cf. 
equation we also see that, together with equation (jA31|) . 

m x m 2 \m\ = \z J rW\ \a\ = . (AS2) 

In order to determine 9' one notes, from equations (|4Hjl and ()45j1 . that 

\z + w\ 2 - \x\ 2 = (c 2 - s 2 ) (m 2 2 - mfj . (A33) 
The analogous relation for Mr is 

c' 2 - s' 2 = » P + 9 » 2 -» m » 2 = ^^^W 

Ml -Ml m\m\{Ml-Ml\ 1 ' 
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